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QUASAR MICROLENSING: WHEN COMPACT MASSES MIMIC SMOOTH MATTER 
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ABSTRACT 

The magnification induced by gravitational microlensing is sensitive to the size of a source relative to 
the Einstein radius, the natural microlensing scale length. This paper investigates the effect of source 
size in the case where the microlensing masses are distributed with a bimodal mass function, with 
solar mass stars representing the normal stellar masses, and smaller masses (down to 8.5 x lO^^Mo) 
representing a dark matter component. It is found that there exists a critical regime where the 
dark matter is initially seen as individual compact masses, but with an increasing source size the 
compact dark matter acts as a smooth mass component. This study reveals that interpretation of 
microlensing light curves, especially claims of small mass dark matter lenses embedded in an overall 
stellar population, must consider the important influence of the size of the source. 
Subject headings: gravitational lensing - microlensing - dark matter halos 



1. INTRODUCTION 

As the light from a distant source shines through 
a foreground galaxy, the gravitational lensing effect 
of the individual stars within the galaxy can in- 
duce rapid flu ctuations in the apparent brightness 
of the source l|Paczvnskil 119861) . Since the detec- 
tion of the first lensing induced variabili ty, seen in 
the quadruply imaged quasar Q2237-f 0305 (|| rwin et aLl 
[1989: Corrigan ct al, 1991) , there has been intense 
theoretical study of this gravitational microlensing, 
revealing that the degree of brightness fluctuations 
is strongly de pendent upon the s ize of the source 
being lensed fKavser et al." '1986'; Wambsgans s et all 
IT99O: Wambsg anss fc Paczynski 1991 ) and the mass 
function of the lensing objects ijLewis fc IrwinI Il996t 
iWvithe fc Turne J 120011) . potentially providing a useful 
probe of both. Further studies have scrutinized the in- 
fiuence of source size and shape fRcfsdal fc St abeTilll99lL 
[1993, 1997; Agol 2003; Mortonson ct al. 2005j, and more 
recent prograrn s with dedicated monitorin g programs 
llWoznia,k et a,].ll200?i (Al^de et a,l. ll200liTI]lan et a,].l 
12003(1 have finally resulted in high quality microlensing 
light curves that can be compared directly with theoret- 
ical models ( Wvithc et al. 2000; Kochanck 2001). 

More recently, gravitational microlensing has been 
proposed as an explanation of the anomalous flux ra- 
tios seen in multiply imaged quasars, with the addi- 
tion of smoothly distributed dark matter broadening 
the wid th of t he magni fication probability distribution 
l|Schechter fc Wambsganss! l200^ . This result appears 
somewhat surprising given a conjecture that has existed 
for a number of years in microlensing studies, namely 
that the shape of a magnification probability distribu- 
tion should be independ ent of the mass funct ion of the 
micro lensing masses ( Wambsganssi [199^ iLewis fc IrwinI 
11996(1 : the smooth matter limit is approached as the mass 
of the individual objects is decreased, but the surface 
mass density remains the same. Clearly, the significant 
differences seen in the magnification probability distribu- 



tions with smooth and compact matter violates this con- 
j ectur e, a position recently confirmed bv lSchechter et al] 
ll20f)l . 

The fall of the microlensing conjecture therefore im- 
plies that there is a fundamental difference in the action 
of compact masses as opposed to smoothly distributed 
mass. But can compact masses appear as smooth mat- 
ter to a source? It is well known that increasing the 
size of a microlensing source smooths out the variability 
of the light curve, but does this imply that to a large 
source, a population of small masses appears as smooth 
matter? While this situation was not examined in the 
work of iSchechter et all l|2004() . it existence means that 
the breakdown of the microlensing conjecture is further 
blurred by the size of the source under consideration. 
Studying this regime, therefore, is the goal of this paper. 
In Section [21 the numerical approach is outlined, while 
Section[2|discusses the influence of a range of source sizes 
and compares the resulting properties of the magnifica- 
tion maps. In Sections 0) and [S] the implications of this 
study and the conclusions are respectively presented. 

2. APPROACH 

In studying the infiuence of gravitational microlensing, 
this study employs the backwards ray tracing algorithm 
of Kavscr et al. (1986) and Wambsganss (1992), whereby 
a large number of rays are fired through a field of mi- 
crolensing masses and collected to form a magnification 
map over the source plane. The important length scale 
for gravitational microlensing is the Einstein radius in 
the source plane, given by 
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where Dij are angular diameter distances between the 
observer (o), lens (l) and source (s) respectively and 
typically this sca le is Vn ^ 0.04pc for cosmological cir- 
cumstances (e.g. i Lewis fc IbatallT998() . As well as the 
mass function of the microlensing bodies, other parame- 
ters need to be specified, including the normalized surface 
mass density (tr) which may be composed of compact and 
smooth masses, and a shear term (7) which describes the 
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Fig. 1. — Examples of the magnifications maps utilized in this paper. Each map is 12 Einstein radii on a side with a normalized surface 
mass density of tr = 0.2 and global shear of 7 = 0.5. The lower-right panel presents the case where half of the surface mass density is in the 
form of solar mass objects and the rest in the form of smoothly distributed matter. In the additional panels, this smooth matter component 
has been replaced with compact mass objects with solar and then sub-solar masses; the individual masses of the sub-solar population are 
noted in the lower left-hand corner of each panel. Note that in each panel, the solar mass stars retain their positions in the lensing plane. 



global asymmetric gravitational lensing influence of the 
overall galactic/intracluster mass distribution. Given the 
potential parameter space, a comprehensive study of po- 
tential microlensing configurations is beyond the scope 
of this paper. As a first step, a representative model was 
chosen with a total surface mass density of cr = 0.2 and 
global shear of of 7 = 0.5. Each simulation covers an area 
of 12^ Einstein radii for a solar mass star and consists of 
2048^ pixels. 

A fiducial model was chosen in which half the surface 
mass density is in the form of solar mass stars, while the 
remaining mass is in the form of smoothly distributed 
matter. In addition to this, 43 additional simulations 
were undertaken. Again, half the mass was in the form 
of solar mass stars, while the remainder is also in the 
form of compact objects, but with progressively smaller 
masses, with each subsequent dark matter mass being 
80% of the proceeding value. The limit of this procedure 
was compact dark matter masses of 8.5 x lO^^M©, re- 
sulting in a total number of microlensing masses in the 
simulation exceeding 4.2 x 10^; it was found that this 
represented the memory limit of available computational 
hardware, while greatly increasing the calculation time 
of the magnification maps, and hence no smaller masses 
were considered. It is important to note that the loca- 
tions of the solar mass stars are identical in all simula- 
tions. 

Figure Q presents examples of the results of this ap- 
proach, with the upper-left panel representing the mag- 
nification maps for the case where half the matter is in 
the form of smooth matter. For the additional panels 
in this figure, this smooth matter component has been 
replaced with compact objects, beginning with IMq and 
continuing down to 1.1 x lO^'^M©. There is a signifi- 



cant difference between the caustic structures apparent 
in the magnification map for the smooth matter case and 
that for all the lensing bodies being of IMq. However, 
as the mass of the compact objects replacing the smooth 
matter component is reduced, there is a very apparent 
evolution in the form of the caustic structure; maps ap- 
pear to possess caustics structure on two distinct scales, 
corresponding to the two microlensing masses under con- 
sideration. For the magnification map with the smallest 
masses, the large scale caustic structure becomes very 
similar to that of the case with smooth matter. However, 
it is also apparent that there is significant small scale 
structure distributed throughout the map. As will be 
seen in the next section, while the magnification maps in 
smooth matter and small mass cases appear similar, the 
statistical properties, especially the magnification prob- 
ability distributions, differ markedly. 

3. ANALYSIS 

To determine the microlensing properties of a partic- 
ular extended source, the magnification maps must be 
convolved with the appropriate surface brightness dis- 
tribution. For the purposes of this study, the source is 
assumed to possess a Gaussian profile with a character- 
istic radius. As the source size is increased, the resulting 
magnification map and magnification probability distri- 
bution are compared with those of the solar mass plus 
smooth matter component. 

Figure El illustrates this procedure for the case where 
the compact dark matter component have individual 
masses of 8.5 x 1O~^M0, the smallest masses that are 
computationally viable at this time. The upper pan- 
els present the unconvolved magnification maps for the 
smooth dark matter (left-hand) and compact dark mat- 
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Fig. 2. — An example of the procedure undertaken in this paper. The left-most column contains a magnification map formed when half 
the surface mass density is in solar mass objects, with the rest smoothly distributed. The centre panel presents the same circumstance, 
except the smooth matter component has been replaced with compact objects each with a mass of 8.5 X IO-'^Mq. Following these are 
light curves for the smooth matter (red) and compact mass (green) cases, taken across the centre of the maps (indicated by the lines 
on the magnification maps). Note that the light curves are in magnitude and the range on the y-axis from +3 to -5 magnitudes about 
the theoretically expected mean. This scale has been omitted for clarity. The right hand column presents the magnification probability 
distributions. The top maps represent a single pixel source, while moving downwards the source size is increased. As it does, the small 
scale caustics structure induced by the smaller masses is washed out and the maps, light curves and magnification distributions become 
the same. While the top panel is unconvolved, the scale-radius of the Gaussian kernel in the subsequent panels are 5.9 X IQ-^, 1.8 X 10^2^ 
5.7 X 10~^ and 1.8 X 10~^ Einstein radii respectively. 



ter (right-hand)^. The maps possess snuilar large scale 
caustic structure, but clearly the presence of small masses 
instead of smooth dark matter has introduced significant 
substructure into the caustic network. Over-plotted on 
these magnification maps are the resultant light curves of 

^ It should be remembered that these maps, by virtue of the 
pixelated nature, are implicitly convolved with a source on the 
scale of a single pixel. 



a source passing horizontally across the middle of the two 
maps and these present quite dramatic behaviour. This 
is reflected in the magnification probability distributions 
which are markedly different. Progressing down the page 
these magnification maps are convolved with larger and 
larger sources. As the source size increases, the sub- 
structure in the magnification map with compact dark 
matter begins to wash out. This is further reflected in 
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Fig. 3. — The solid line in this figure denotes the Einstein radius 
for the smaller mass microlensing simulations, whereas the points 
present the radius of the Gaussian kernel at which the difference 
in the binned, cumulative magnification patterns differed by less 
than / = 0.05 (circles) and / = 0.10 (squares); note the relation 
shows an linear trend in this logarithmic space. The discrepant 
points are due to "noise" in the binned distribution preventing 
the fit from reaching the fitting criteria. The arrows in this figure 
denote the magnification maps which did not achieve this criterion 
and the magnification distributions are too different to allow an 
equivalence between the two maps even with smoothing on the 
scales considered. 



the subsequent light curves, with the degree of variabihty 
decreasing. Interestingly, as the source size is increased, 
the light curves for the compact dark matter and corre- 
spondingly convolved smooth matter cases become very 
similar, again a fact seen in the magnification probabil- 
ity distributions which too become identical. With even 
larger smoothing, the magnification maps remain identi- 
cal and hence the compact matter will effectively appear 
as a smooth dark matter component. 

To examine this further, the entire sample was sub- 
jected to the same convolution procedure, starting with 
small Gaussian sources and increasing the radius of the 
kernel until the magnification probability distributions 
became similar. In defining this threshold, the binned 
cumulative magnification distributions were compared, 
with the difference critical convolution radius defined as 
being that at which 
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where and are the cumulative probability distri- 
butions; such cumulative distributions ease differences 
in the differential magnification probability distributions 
due to noise. The factor / represents the fractional dif- 
ference between the two magnification distribution and 
for the purposes of this study we choose / = 0.05 and 
/ = 0.10 as representative values. Of course, the choice 
of the value of / is rather arbitrary and other choices 
could be made. Hence, the functional form employed 
here should be viewed as illustrative rather than defini- 
tive. The results of this procedure are presented in Fig- 
ure 121 here the solid line represents the Einstein radius 
of the smaller mass component, with the points denoting 



the critical radius as defined by the above criterion. This 
figure shows that there is a clear trend in this statistic 
with the critical Gaussian radius increasing in step with 
the Einstein radius. For / = 0.05 and mass scales from 
^ IO^^'Mq up to '^O.1M0, the ratio of these two quan- 
tities is (remarkably) constant at ~ 4 (with a couple of 
discrepant points due to noise spikes in the cumulative 
distribution delaying convergence), although this ratio 
drops to ~ 1.1 for / = 0.10. Note, the arrows denote 
masses at which smoothing on the scale of less than 5 
Einstein radii failed to result in convergence between the 
two magnification distributions. A reexamination of Fig- 
ure ^ illustrates that at this high mass end, the caustic 
structure introduced by the compact masses is similar in 
scale to the overall caustic network and it is understand- 
able that this structure has not been smoothed out on 
the scales under consideration. 

Returning to Figure |3 it is clear that the third and 
fourth row of panels straddle this critical region, with 
magnification distributions for the smaller Gaussian 
source appearing quite different for the smooth matter 
and compact cases. Considering the lower panel, and 
hence larger source, the distributions in both scenarios 
have become very similar. Notice in this case, even when 
smoothed with the larger source, the resulting magni- 
fication distributions possess significant structure which 
is refiected in quite dramatic events in the microlensing 
light curve. Hence the origin of the similarity in the 
distributions is not due to all of the magnification struc- 
ture being smoothed out of the microlensing maps. For 
cases where the difference in mass between the compact 
dark matter and the solar mass stars is less, then a lot 
of the structure in the magnification maps is smoothed 
out, with the resulting light curves possessing very little 
structure. 

3.1. Further simulations 

So far this paper has considered a single combination 
of microlensing optical depth and shear, but how gen- 
eral are the results uncovered thus far? In addressing 
this question, three more sets of microlensing parame- 
ters were employed, namely ((7,7) = (0.2,0.0), (cr, 7) = 
(0.2,0.2) and ((7,7) = (0.6,0.6). As with the previous 
simulations, the source region under consideration cov- 
ered 12^ Einstein radii for a solar mass star and consisted 
of 2048^ pixels. Again, in each case, half of the optical 
depth comprised of solar mass stars while the other half 
was in the form of smooth matter or compact objects, of 
mass either 3.5 x IO^^Mq or 1.2 x lO'^Mg. 

Rather than simply repeating the analysis presented in 
Section |3 a slightly different approach was undertaken. 
As revealed in Fig. as the masses of the compact ob- 
jects decreases, the structure in the magnification map 
tends to that of the smooth matter case, with the smaller 
masses providing strong, but localized, perturbations, 
and it is these perturbations which are washed out by 
convolving with a large enough source; note, however, if 
the masses of the compact objects are large enough, they 
produce gross changes into the magnification map such 
that convolving cannot make the magnification identi- 
cal to the smooth matter case. Hence, given the small 
masses considered in this second set of simulations, we 
would expect that the structures in the compact matter 
and smooth matter magnification maps should become 
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Fig. 4. — The variation in the su btra cted convolved magnifica- 
tion maps, as discussed in Section 13.11 The differing line styles 
correspond to the three parameter sets employed (as denoted by 
the key) where as the lighter lines represent the compact masses 
of 3.5 X 10~^Mq, whereas the heavier lines are for masses of 
1.2 X 1O~^M0. Note that the cases with a = 0.2 possess simi- 
lar structure and overlap in this figure. 



similar, if convolved with a large enough source. 

The analysis procedure again consisted of convolving 
the smooth matter and compact matter magnification 
maps with a Gaussian profile. A residual map was con- 
structed by converting the smoothed magnification maps 
to magnitudes and then subtracting them. These residu- 
als possess a Gaussian-like profile, centred upon zero, and 
so a residual dispersion, <7Mag was calculated. The re- 
sults of this procedure are presented in Figure^ in which 
this residual dispersion is plotted against the smooth- 
ing radius, for the simulations under consideration. The 
trend is as expected, with the residual map showing large 
variations, with crMag ~ 1 at small smoothing radii for 
((7,7) = (0.6,0.6), but as the smoothing radius is in- 
creased, aMag falls, and hence the convolved magnifi- 
cation maps present similar structure. Also, aMag falls 
faster for the smaller compact masses, as their magnifi- 
cation structure gets washed out at smaller radii. 

Interestingly, the curves for the larger mass compact 
objects in Fig. 01 clearly show a plateau-like structure, 
remaining constant over a range of smoothing radii, and 
then beginning to decrease. Furthermore, the similarity 
of the curves for both masses considered suggests that 
the smaller mass component possesses a similar plateau 
structure, existing in the regime below the minimum 
smoothing radius we considered. Hence, this again re- 
veals that there appears to be a critical radius below 
which the influence of the compact matter changes char- 
acter, effectively representing a smooth mass component. 
A systematic study of this critical radius for a range of 
microlensing parameters, and the influence of physical 
size and structure of the emission regions in quasars, will 
be the subject of a forthcoming contribution. 



4. IMPLICATIONS 

One immediate implication of the results presented in 
this paper is that if the dark matter is in the form of 
relatively small compact masses, then the resulting sta- 
tistical properties of a microlensing light curve will be de- 
pendent upon the size of the source under consideration. 
It is important to note that this is somewhat different 
to the usual size dependence of microlensing statistics 
where there is a straightforward smoothing of the light 
curve as the source size is increased, with a resultant 
narrowing of the magnification probability distribution. 
Rather, a small source and large source, where relative 
sizes of small and large can be determined from the rela- 
tionship presented in Figure |3| will be subject to signif- 
icantly different magnification patterns and correspond- 
ing statistics. 

The effect influences the question of how well we can re- 
veal the nature of the dark matter component, if present. 
In realistic gravitational microlensing scenarios, the dis- 
tribution of matter in the lensing galaxy will not be 
a clean two-mass component population, but rather a 
certain mass function (or even a combination of mass 
functions). Depending upon the source size for a given 
system, this mass function will suffer a cutoff point be- 
low which observations will not be able to distinguish 
between smooth and granular matter. This does not 
mean that small lens masses and large sources combi- 
nation with a single mass component cannot produce 
microlensing imprints. Similarly means that for large 
enough source sizes, sub- Jupiter masses cannot be de- 
tected unambiguously and that claims involving plane- 
tary and sub-planetary populations that produce light 
curves fluctuations should be thoroughly justified (i.e. 
fCoUev & Schild 2003). For example, from Fig.^l for / = 
0.05, we see that for microlenses masses of ~0.04 Mq, 
there are still source sizes (~0.7 ER) that make such a 
low-mass population to appear as a smooth distributed 
matter component. Considering the corresponding phys- 
ical lengt h scale in microlense d quasars corresponds to 
~ 0.04pc l)Lewis fc IbatallT998l) . it is too large to repre- 
sent the optical/UV emitting region of quasars. However, 
for more realistic quasar accretion disk sizes of <0.1 ER, 
"secondary" populations with masses <10~^ Mq will be 
completely smoothed out. Since the trend in Fig. O is 
quite constant for both considered values of /, this ap- 
proach would also be useful in putting limits to source 
sizes, though a better understanding of the role of the 
mass populations is needed. 

To some extent, it seems that our results are in contra- 
diction with those obtained by Schechter et al. (2004). In 
that work, the authors concluded that "the magnification 
probability distribution for two disparate components is 
not exactly that of a single component and a smooth com- 
ponent" . The apparent contradiction vanishes when it is 
realised that Schechter et al. (2004) did not account for 
the source size effect (or, in other words, they assumed 
the same pixel size source in their study). Indeed, this 
is an important issue, because observationally the source 
effect will be always present. This means that the ad- 
ditional structure in the magnification maps induced by 
the low-mass component, seen by Schechter et al. (2004), 
is again blurred out by the source. So, the dependence 
of the magnification probability on the higher-order mo- 
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ments of the mass distributions is no longer true when 
introducing the source effect (at least for certain regions 
of the mass distribution function) . 

5. CONCLUSIONS 

This paper has presented a study of the influence of 
source size on the properties of gravitational microlens- 
ing with a bimodal distribution of lensing masses. Given 
the computational limitations, the study focused upon 
a set of macrolensing parameters, but considered the 
influence of compact dark matter on a range of mass 
scales. It was found that small sources resolved the caus- 
tic network produced by the smaller masses, but as the 
source size was increased this fine scale network becomes 
washed out. When compared to the scenario where the 
smaller masses are replaced with an equivalent quantity 
of smooth matter, it is seen that the light curves of the 
smallest sources in both circumstances are quite differ- 
ent. As the source size is increased, the light curves, and 
corresponding magnification probability distributions, in 
each case become identical, indicating that there is a crit- 
ical source scale where compact dark matter behaves, in 
a gravitational lensing sense, like smoothly distributed 
matter. This critical source scale is dependent upon the 
ratio of the stellar to dark matter masses, with a rel- 
atively linear trend in log-log space (i.e. source sizes 
smaller than the circles in Figure |21 are small enough 
to resolve the compact dark matter, and we can expect 
to see their signature in the microlensing light curve, 
whereas for source size greater than the circles, this sig- 
nature will be smoothed out). 

The results found in this paper point to another prob- 
lem concerning magnification probabilities. To what ex- 
tent magnification distributions are showing valuable in- 
formation in real gravitational lensed systems? In other 
words, can we characterize systems according to their 
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magnification distributions? To answer these questions 
we would need to explore a large range of combinations of 
a and 7, together with different mass distribution func- 
tions; while computationally expensive, we are currently 
planning the first stages of such an exploration. How- 
ever, it seems clear that when the source effect is taken 
into account, the information available from the mag- 
nification distributions can be rather poor. Of course, 
the magnification distributions are only the zeroth level 
statistic that can be used to examine the influence of 
gravitational microlensing, but clearly the results pre- 
sented in this contribution similarly influence higher or- 
der temporal statistics that are applied to microlensing 
light curves. 

In closing, this contribution returns to the question of 
the microlensing conjecture that was discussed earlier. 
The fall of the conjecture implies that the mass func- 
tion of the compact objects was potentially amenable 
to microlensing observations as it is now apparent that 
this will imprint a signature into any microlensing light 
curves. However, the issue has been clouded further with 
the results presented in this paper illustrating that the 
size of a source being of vital importance when address- 
ing the question of how well the properties of an underly- 
ing mass function can be determined. While introducing 
a mass function expands the potential parameter space 
enormously, it is important to understand the influence 
of the source on the values of microlensing statistics and 
hence future contributions will investigate its role in more 
general microlensing scenarios. 
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